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We analyze the applicability of the Rosenfeld entropy scaling relations to the systems with the 
core-softened potentials demonstrating the water-like anomalies. It is shown that the validity of the 
of Rosenfeld scaling relation for the diffusion coefficient depends on the thermodynamic path which 
is used for the calculations of the kinetic coefficients and the excess entropy. In particular, it is valid 
along isochors, but it breaks down along isotherms. 
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In 1977 Rosenfeld proposed the relations connecting 
transport properties of a liquid with the excess entropy 
[l[ . In order to write down these relations one should use 
the reduced forms of the transport coefficients: 



D* = D 



V =11- 
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where D and 77 are the diffusion coefficient and the viscos- 
ity. According to the Rosenfeld suggestion, the reduced 
transport coefficients can be expressed in the form 
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where Sex = {S — Sid)/ [NkB) is excess entropy of the 
liquid, X is the transport coefficient, and ax and bx are 
the constants which depend on the studying property Q ■ 
Interestingly, the coefficients a and h show an extremely 
weak dependence on the material and can be considered 
as universal. 

In his original works Rosenfeld considered hard 
spheres, soft spheres, Lennard- Jones system and one- 
component plasma [H, 0] . After that the excess entropy 
scaling was applied to many different systems including 
core-softened liq uids liquid metals [1, 0], binary 

mixtures 11 1, ionic liquids 12, 13], network- forming 
liquids fsl. Il2l| water [3], chain fluids [lH and bounded 
potentials [TL iH 17|. 

Nevertheless, controversies still remain. For example, 
up to the moment it is not clear whether the Rosenfeld 
scaling relations are applicable to the core softened sys- 
tems. Some publications state that the scaling relations 
are valid for such systems [H, 0] , while in our recent work 
it was shown that the scaling relations may break down 
for the core softened systems Q- This article presents a 
discussion of this contradiction. Basing on the molecular 
dynamics simulations of the two core-softened systems 
we show that the validity of the Rosenfeld scaling rela- 
tions depends on the thermodynamic path which is used 



for the calculations of the kinetic coefficients and the ex- 
cess entropy. In particular, the exponential functional 
form of the relation between the excess entropy and the 
reduced diffusion coefficient holds along isochors while 
along isotherms one observes its breakdown. 

Two systems are studied in the present work. The 
first one is the core-softened system introduced by de 
Oliveira et al [3] ■ This system is described by the spher- 
ically symmetric potential represented by a sum of a 
Lennard- Jones contribution and a Gaussian-core inter- 
action (LJG): 



U{r) = 4e 



12 



- ae ■ exp 



1 



(4) 

with a ~ 5.0, tq/u = 0.7 and c = 1.0. This model can 
qualitatively reproduce water's density, diffusivity, and 
structural anomalies. The diffusivity of this system was 
studied, for example, in the papers [y, |l8|. In the papers 
0, @] the validity of the Rosenfeld scaling relation for 
this system was checked. In this paper we discuss these 
results and compare them with our calculations. 

The second system studied in this work is a repulsive 
shoulder system (RSS) introduced in our previous work 
[1^. This system has a potential 



U{r) = (^)" + ie- [l-tanh(A:o{r-fTi})], (5) 

where a is the "hard" -core diameter, cti ~ 1.35 is the 
soft-core diameter, and k^ = 10.0. In Ref. [20] it was 
shown that this system demonstrates anomalous ther- 
modynamic behavior. In our previous publication Q the 
Rosenfeld relation for this system was studied. It was 
shown that the scaling relation for the diffusion coeffi- 
cient breaks down for this system in the anomalous dif- 
fusion region. 

In this paper we use the dimensionless quantities: f = 
r/cr, P = Pa^/e, V = V/Na^ = \/p,f ^ ksT/e. Since 
we use only these reduced units we omit the tilde marks. 

For the investigation of the LJG potential we simu- 
late a system of 1000 particles in a cubic box for the 
densities ranging from p — 0.01 till p = 0.35 with the 
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FIG. 1: (a) The diffusion coefficient of the LJG system along 
some isotherms, (b) The excess entropy of the LJG system 
for a set of isotherms. (Color online) 



Step 5p = 0.01. The time step used is dt = 0.001. The 
equilibration period consists of 1 • 10^ time steps and the 
production period - 2.5 • 10^ time steps. During the equi- 
libration the temperature is kept constant by velocity 
rescaling while during the production cycle NVE-MT) is 
used. The equations of motion are integrated by velocity- 
Verlet algorithm. The following isotherms are simulated: 
T = 0.2; 0.3; 0.4; 0.5; 0.6; 1.0 and 1.5. 

The simulation setup of the RSS was described in Ref. 







The excess entropy in both cases was computed via 
thermodynamic integration method. For doing this we 
calculate excess free energy by integrating the equa- 
tion of state along an isotherm: J^^t ~ ^Nk^T ^ 
_2_ jp p{p )-p ksT rpj^g excess entropy is computed 

via Sex = ^NksT ■ 



FIG. 2: The densities of maximum and minimum of the dif- 
fusivity and the excess entropy for the several isotherms. 



Fig. [Ija) presents the diffusion coefficient of the LJG 
system along a set of isotherms. One can see that at low 
temperatures the diffusion coefficient shows nonmono- 
tonic behavior which is known as the diffusion anomaly. 
Fig. [ijb) demonstrates the excess entropy for the same 
set of isotherms. From these figures it follows that 
the diffusion coefficient and the excess entropy along an 
isotherm have similar qualitative behavior. However, the 
location of the extremum points of diffusivity and ex- 
cess entropy is different (Fig. [2])- It means that there are 
some regions where one function increases while another 
one decreases and vice versa. Clearly, this kind of be- 
havior can not be consistent with the Roscnfeld scaling 
formula. From this one can conclude that the Rosenfeld 
scaling relation is not applicable to the diffusivity along 
an isotherm for the LJG model. 

Fig. 13] presents the logarithm of reduced diffusion coef- 
ficient vs excess entropy along a set of isotherms. One can 
see that at low temperatures the curves demonstrate a 
selfcrossing loop like the one observed for the RSS model 
in Ref. [7|. This loop becomes less pronounced with 
increasing the temperature and at T = 0.5 the self cross- 
ing disappears. At high temperatures the curve comes 
to the straight line limit which corresponds to the small 
influence of the soft core at such temperatures. 

It can be seen from Figs, [ija) and (b), that both 
the diffusivity and the excess entropy along isochors 
are monotonous (this corresponds to the vertical lines 
through the data points). It allows to expect that the 
exponential relation between the diffusion coefficient and 
the excess entropy holds along isochors. Fig. |4] shows the 
ln{D*) vs Sex along a set of isochors. From this figure 
one can see that all curves with good accuracy corre- 
spond to the Rosenfeld scalingrelation. This result is in 
agreement with with Refs. [1,1^. 

In order to show that the dependence of the Rosenfeld 
scaling relations on the thermodynamic path is a general 
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FIG. 3: The logarithm of the reduced diffusion coefficient 
for the LJG model along the several isotherms. The insets 
correspond to the temperatures T — 0.3 (upper panel) and 
T = 0.5 (lower panel) in the larger scale. (Color online) 
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FIG. 4: The logarithm of the reduced diffusion coefficient for 
the LJG model along the several isochors. (Color online) 



property of the core-softened systems we study another 
core softened potential introduced above - the repulsive 
shoulder system (RSS). In our previous publications we 
showed that this system demonstrates anomalous dif- 
fusion behavior at low temperatures (20l | and that the 
Rosenfeld scaling relation along isotherms is not applica- 
ble Q : the curves demonstrate the self crossing loops like 
in the case of the LJG system. Fig.[5]shows the logarithm 
of the reduced diffusion coefhcient along a set of isochors 
for the RSS. As it can be seen from this figure, the de- 
pendence of D* on Sex is linear. However, the slope of 
the line shows an isochor dependence. Fig. [S] shows that 



FIG. 5: Reduced diffusion logarithm for RSS along a set of 
isochors. 

the slope remains approximately constant for low densi- 
ties (p = 0.3 — 0.55) while on increasing the density the 
slope also increases. 

The diffusion coefficient along isotherms as a function 
of the density and the reduced diffusion coefficient as a 
function of excess entropy are shown in Figs.[6l[a) and (b). 
As it can be seen from Fig. IDJa), the density p = 0.55 
corresponds to the density maximum at T = 0.2. In 
Ref. [19] it was shown that this system is effectively 
quasibinary. The reason for the complex behavior of this 
system is related to the competition between two length 
scales a and cti. This allows to suggest that the density 
p = 0.55 belongs to the boundary between two different 
regimes. It seems that the densities from p = 0.6 till 
p — 0.8 belong to a cross section region between low- 
and high-density regimes which leads to the shift of all 
curves with respect to each other. One can expect that at 
higher densities the logarithm of D* comes to a straight 
line again but the slope of this line should be different 
from that at low densities. 

In conclusion, in the present article we carry out a 
study of applicability of the Rosenfeld excess entropy 
scaling relation for the diffusion coefficient in the case 
of the systems with the core-softened potentials. We 
find that the validity of the of Rosenfeld scaling rela- 
tion depends on the thermodynamic path. In particular, 
it is valid along isochors (although the scaling parame- 
ters show the isochor dependence), but it breaks down 
along isotherms. Interestingly, both the excess entropy 
and the diffusivity are monotonous along isochors while 
along isotherms they show the anomalous increase at low 
temperatures. It is this inconsistency in the location of 
the minimum and maximum points of the diffusion coef- 
ficient and the excess entropy along the isotherms which 
leads to the breakdown of Rosenfeld scaling relation for 
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FIG. 6: (a) The diffusion coefScient for the RSS system along 
several isotherms. The inset enlarges the isotherms T = 0.2 
(squares) and T — 0.25 (circles), (b) The logarithm of the 
reduced diffusion coefficient as a function of Sex along a set 
of isotherms for RSS. (Color online). 



these systems. This allows to suggest that the scaling 
relation is valid for the core-softened systems along any 
thermodynamic path if both Sex and D* are monotonous 
along it. However, this suggestion requires further inves- 
tigation. 

Another class of systems which attract a lot of at- 
tention with respect to the Rosenfeld scaling relation is 

bounded potential systems 0, Hm. It is important 

to note that the conclusion of this paper is not appli- 
cable to these systems. For example, in Refs. [l6|, |l7| 
scaling relation for the Gauss Core Model ([2l[) along 
isochors was studied. According to these publications, 
the curves deviate from exponential dependence. This 



allows to suggest that the relation between the dynamic 
and thermodynamic properties for the bounded potential 
systems is different from the systems with the hard core. 
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